Spatial pair correlations of atoms in molecular dissociation 
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We perform first-principles quantum simulations of dissociation of trapped, spatially inhomoge- 
neous Bose-Einstein condensates of molecular dimers. Specifically, we study spatial pair correlations 
of atoms produced in dissociation after time of flight. We find that the observable correlations may 
significantly degrade in systems with spatial inhomogeneity compared to the predictions of idealized 
uniform models. We show how binning of the signal can enhance the detectable correlations and 
lead to the violation of the classical Cauchy-Schwartz inequality and relative number squeezing. 
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Producing and utilizing quantum mechanical correla- 
tions between entangled particle pairs are major themes 
in quantum mechanics research. In this paper we investi- 
gate correlations of the kind used by Einstein, Podolsky 
and Rosen (EPR) to argue that local realism is inconsis- 
tent with quantum mechanics being a complete theory 
of nature [l| . EPR type non-local correlations were first 
produced experimentally in the form of photon pairs from 
nondegenerate parametric down-conversion [2]. Here, the 
correlations are observed between the optical quadrature 
amplitudes of spatially separated signal and idler beams, 
which are analogous to canonical position and momen- 
tum variables originally discussed by EPR. 

It should be possible to extend this type of experiment 
to ultra-cold atoms, which would allow the investiga- 
tion of EPR correlations for particles with non-zero rest 
mass. These might be produced by the transfer of cor- 
relations from photons to atoms or else by nonlinear 
matter- wave interactions. The latter include atomic four- 
wave mixing via Bose-Einstein condensate (BEC) colli- 
sions 0, Q and matter- wave down-conversion via dissoci- 
ation of a molecular BEC 0, 0, B Hi • Although the direct 
measurement of matter-wave quadrature correlations of 
this type is a challenging task [9| , a simpler measurement 
of density-density correlations - as a prerequisite for the 
EPR quadrature correlations - can be performed with 
current experimental techniques. 

In this paper we perform first-principles quantum sim- 
ulations of dissociation of a BEC of molecular dimers into 
pair-correlated bosonic atoms. This has been demon- 
strated experimentally using ultracold (but not Bose 
condensed) ^''Rb2 dimers [ifl], although the atom cor- 
relations have not yet been measured. In the case of 
fermionic atoms 11, l3|, such as in dissociation of 
dimers [l3| , the density-density correlations between the 
atoms in two different spin states have been measured; 
this case, however, currently resists first-principles quan- 
tum simulations. Nevertheless, many insights from the 
present bosonic case also apply to the fermionic case. 



and to correlations produced via BEC four-wave mixing 
fl3|. For example; the importance of operationally well 
defined measures of correlations, the role of spatial in- 
homogeneity, and the trade-off between image resolution 
and binning to enhance correlation detection. 

Measurements now challenge theory to give precise 
quantitative descriptions of the observable correlations. 
Experiments measuring atom shot noise in absorption 



images II 
tectors 11 




and using microchannel plate de- 



^ ^ have demonstrated atom-atom correlations 

1^ . 14 1 . However, their quantitative theoretical analysis 



is lacking. There is also a need to determine the best op- 
erational measures 12|, |l8| of correlations for these exper- 
iments. The present work is a step towards these goals, 
and a benchmark for other (approximate) theories. 

Previous theoretical work on molecular dissociation 
has explored the atomic correlations in simplified cases 
such as ID systems 0, without depletion of the 
molecules during conversion [l^ , or for spatially uniform 
systems 
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18j . The present work builds on the charac- 
terization of the atomic correlations in momentum space 
in Ref. [lH. Here, we extend the analysis to spatial cor- 
relations and study molecular dissociation in 3D under 
experimentally realistic conditions, including the effects 
of inhomogeneity of the initial molecular BEC [3] and 
time-of-flight expansion. We focus on the use of absorp- 
tion images, obtained by line of sight integration, to mea- 
sure the atomic correlations. 

To investigate the quantum dynamics of molecular dis- 
sociation we use stochastic equations in the positive P- 
representation 20). In this method, one simulates the 
evolution of four complex stochastic (c-number) fields 
\I'i(x, i) and <i>i(x, f) [<l>*(x, t) ^ *i(x, t)], representing 
the field operators ^i(x, t) and ^'|(x, t), respectively, 
where i = stands for molecules and i = 1 for atoms. 
Averages of the stochastic field products over a large 
number of trajectories (. . .)st correspond to quantum 
mechanical ensemble averages of normally-ordered oper- 
ator moments. For example ([>I'|(x, t)]'^[>I'j(x', t)]") = 
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([$,(x,t)]'=[*,(x',t)]")3t. 

The stochastic differential equations governing the 
quantum dynamics of dissociation are [18 1 



responding atom number operator is given by 
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Here, (j{x,t) {j — 1,2) are real, independent Gaus- 
sian noises with zero means and nonzero correlations 
(0(x,t)a(x',i')).t = S,kS{x - x')dit - t'), mi and 
mo = 2mi are the atomic and molecular masses, A 
is the detuning corresponding to the energy mismatch 
2?iA = 2i?i — Eq between the free two-atom state at 
the dissociation threshold and the bound molecular state, 
and X is the coupling responsible for coherent conversion 
of molecules into atom pairs, e.g., via Raman transitions 
or a Feshbach resonance (see Refs. 



18|,l21|] for details). 



We have ignored s-wave scattering interactions because 
they do not significantly affect the dynamics or the cor- 
relations for durations resulting in less than ~ 10% con- 
version [1^. The characteristic timescale for this can be 
estimated using a Fermi's golden rule calculation 
giving t < 0.ir-\ where F = x'^imi/2h)^/'^ y^l/n is 
the molecular decay rate. The duration of dissociation 
should also be shorter than the timescale for losses due 
to inelastic collisions [2^. 

We assume that the molecular BEG is initially in a co- 
herent state, whereas the atoms are in the vacuum state. 
Once the dissociation is suddenly switched on [2^, the 
trapping potentials are simultaneously switched off, so 
that the evolution takes place in free space. For molecules 
at rest, the excess of potential energy is converted into ki- 
netic energy, 2?i|A| 2?i^fc^/(2mi), of dissociated atom 
pairs with equal but opposite momenta around ikp, 
where |ko| = •y/2mi | A|/fi,. This is the physical origin of 
the expected correlations between the atoms. Ideally, in 
the time-of-flight expansion the momentum correlations 
are converted into position correlations between diamet- 
rically opposite atoms in the far field. In non-ideal cases, 
such as in strongly inhomogeneous systems with large 
momentum uncertainty, or with insufficient time for ex- 
pansion, the spatial correlations may degrade. 

In absorption imaging the number of photons detected 
by each camera pixel determines the number of atoms 
contained in the volume of a narrow column in the imag- 
ing laser's propagation direction, which we denote as z. 
Denoting the area of the camera pixel about the point 
r = x_L = (x,y) on the detection plane as ^(r), the cor- 
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On a computational grid, the number operator is 
related to the integrated 2D column density n±(r) = 
/ dzn(x) via = n±{r)AxAy, where n(x) — 

^i(x)^i(x) is the 3D density. Ax and Ay are the lattice 
spacings, and A(r) = AxAy. 

Gorrelation between the atom number fiuctuations in 
a pair of different pixels can be quantified via the nor- 
malized number-difference variance (r ^ r') 



Vr. 



{[A{N^-N^-)?) 



= 1 



(: [A(jV,~jV,0]^:) 



(3) 

where AC ~ C — (C) is the fluctuation in C and the 
colons :: indicate normally-ordered operator products. 
This definition uses the conventional normalization with 
respect to the shot-noise level of Poissonian statistics, 
such as for a coherent state, (Nj.) + {Nr') . Vr^r' = 1 for un- 
correlated signals. Variance smaller than one, Vr,r' < 1, 
implies reduction (or squeezing) of fluctuations below the 
shot-noise level and is due to correlation between parti- 
cle number fluctuations in the r and r' pixels. Perfect 
(100%) squeezing of the number-difference fluctuations 
corresponds to Vry = 0. 

The number-difference variance is related 
to Glauber's second-order correlation function 
5(2)(r,r') = {:n^{r)n^{v'):)/{{n^(v)){h^ir'))). In 
the simplest symmetric case, with {n±{r)) — {n±{r')) 
and (7*^^-'(r, r) = g(^^(r',r'), the relationship is 



(4) 



In order that T4,r' is suppressed below the level of clas- 
sical, uncorrelated statistics, Vry < 1, it is necessary 
and sufficient to have i?'-^-' (r, r') > ^'^^^(r, r) which is a 
direct signature of the violation of the classical Gauchy- 
Schwartz inequality [25| . Similarly, Vry > 1, implies 
g(^)(^r,r') < f/'^^^(r,r) which agrees with the classical 
Gauchy-Schwartz inequality. Thus, Eq. ([4]) shows that to 
have a non-trivial quantum correlation it is not enough 
to have g'^^\r,r') > 1, but one has to have (7'^'(r,r') > 
g(^^(r, r). Note that a non-ideal atom detection efficiency 
adds a factor of 77 < 1 to the second term in Eq. (|3|) , but 
does not affect the Gauchy-Schwartz inequality. 

An immediate consequence of Eq. ^ is that the ob- 
served number-difference squeezing can be insignificant 
if the pixel occupation number is very small due to a 
small pixel size, even if the pair correlation f/^^-'(r,r') is 
much larger than (7(^^(r,r). In other words, fine spatial 
resolution can degrade the number-difference squeezing 
between the elementary pixels, and larger pixel size or 
binning favor observing strong squeezing. 
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FIG. 1: (color online) Dissociation of a molecular condensate 
in free space, for dissociation time of t = 0.2 ms with x ~ 
7 X 10"'^ m^/^ s"^ and a further expansion of te, — 2.5 ms 
with X = [261] . Each dimension is 55 /im long and contains 
256 lattice points. 10, 000 stochastic paths were averaged, (a) 
Slice through z = of the 3D atomic density n(x) = (n(x)) 
after expansion, (b) 2D column density n^{v) = (n^(r)) af- 
ter integration along the z axis, (c) local pair-correlation, 
g'^^\v,v). The g'^-* data is too noisy outside the dissociation 
sphere where the signal is negligible and should be discarded 
at radii r > 15 /^m. (d) Pair-correlation at opposite locations, 
(;'^'(r, — r). (e) Pair-correlation (;'^'(r, rgo) at orthogonal lo- 
cations, (f) Number-difference variance Vt,-t at opposite lo- 
cations, (g) Binned atomic signal on the detection plane, for 
bins of size 32 x 32 pixels, and (h) number-difference variance 
between the opposite bins. 



Fig. [T] shows the results of a simulation of the dissoci- 
ation of a small molecular BEC of size '^1.4 /im, con- 
taining about 205 molecules (in subsequent examples we 
simulate larger condensates, containing 874 and 9.2 x 10^ 
molecules [26[ ) . The images relate to the atomic field af- 
ter dissociation and expansion. There are about 30 atoms 
in total in this example, corresponding to 7% conversion. 



Fig. [Ija) shows a slice through the 3D atomic density, 
which forms an expanding spherical shell about 14 fim. in 
radius. Fig.[T]Jb) shows the 2D column density after inte- 
gration along the z axis. The projection of the spherical 
shell of atoms is clear and is similar to the experimental 
observations of Refs. [3, [l^]. Figs. [Hc)-(d) show cor- 
relations in column densities at the same, opposite, and 
orthogonal locations: (c) shows the local pair correla- 
tion, f;'^^^(r,r), which has a thermally bunched character 
in 3D, 5^^-'(x, x) = 2 [3], while the obtained value of 
(7'^^(r, r) ~ 1.2 < 2 is due to the integration along z; 
(d) shows strong pair-correlation at opposite locations, 
f;*-^-'(r, — r) ~ 2.05 > .g'^^^(r, r), originating from the mo- 
mentum correlations of dissociated atom pairs; and (e) 
shows an uncorrelated signal of (;^^^(r,rgo) ~ 1, where 
rgo corresponds to a 90° rotation of the original image 
about the origin. In (a)-(f), the elementary pixel size is 
given by the computational grid, with 0.215 /itm spac- 
ing in each dimension. The numerical convergence of our 
simulations is ensured by smaller grid sizes reproducing 
the results within the stochastic sampling errors. 

In the example of Fig. [1] we have a violation of the 
classical Cauchy-Schwartz inequality as g^^^ (r, — r) > 
(j,(2)(j. J.) However, due to small occupation numbers 
of the elementary pixels, (iVr), the associated number- 
difference variance at the diametrally opposite locations, 
Vr,-r, shows very little squeezing in Fig. [^f), as ex- 
plained by Eq. @] it deviates from K.-r = 1 by only 
— 1.7 X 10"'^ on the projected dissociation sphere where 
the signal is maximal. Figs. [l{g) and (h) show that mea- 
surement bins of size 32 x 32 pixels improve the number 
difference squeezing to V^-r — 0.4 (60% squeezing). 

To understand the dependence of the correlation 
strength on the initial size of the molecular BEC and 
the expansion time, we analyzed the binned number- 
difference variance for three different cases. Fig. [2] shows 
the dependence of the angle averaged variance ^rc-ro 
on the expansion time, for different bin sizes. At each 
time, the angle averaged result is at the radius tq of the 
dissociation sphere, where the signal is maximal. 

As we see from Fig. [2l the squeezing improves with 
larger bin sizes in all cases, consistent with the experi- 
mental results of Ref . [l3| . For the case of a small con- 
densate, Fig. HUa), the squeezing degrades with the ex- 
pansion time due to the large momentum uncertainty of 
the initial molecular BEC. The resulting center-of-mass 
momentum offset of correlated atom pairs, causes them 
to fail to appear in diametrically opposite bins in the far 
field. Increasing the bin size captures the pairs in the 
opposite (larger) bins and restores the pair correlations. 
For a larger molecular condensate, as in Fig. [2{c), the 
complimentary effect is that the position uncertainty may 
give the atom pairs a center-of-mass position offset, again 
preventing them from appearing in the diametrically op- 
posite bins in the near field. In this case, the squeez- 
ing improves with expansion, corresponding to a more 
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FIG. 2: (color online) The effect of binning on the angle av- 
eraged number-difference variance between the opposite bins, 
Vro,-TQ, as a function of the expansion time following 0.2 
ms of dissociation. Graphs (a), (b) and (c) correspond, re- 
spectively, to the initial molecular BEG size of 1.4, 3 and 6.4 
fim [26]. The different curves are for different bin sizes. The 
largest bins into just four quadrants on the detection plane. 

complete conversion of the intrinsic opposite-momentum 
correlations into spatial correlations in the far field. The 
case of Fig. [2jb) is intermediate and is affected by the 
competition between these two effects, resulting in the 
optimum expansion time and implying the existence of 
an optimum size of the molecular BEC for a given ex- 
pansion time. 
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